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ABSTRACT 

Context. Theory and observations of heavy element nucleosynthesis are in conflict with one-another. Theory states that in the most 
metal-poor stars, the rapid (r-) neutron-capture nucleosynthetic process would be dominant over the slow (s-) process. The most recent 
determinations of r- and s-process yields do not support this. 

Aims. We provide measurements of the Ba isotopic fractions for five metal-poor stars derived with a local thermodynamic equilibrium 
(LTE) analysis with ID model stellar atmospheres. This increases the comparisons with heavy element nucleosynthesis theory. 
Methods. We use high resolution (R = A/AA = 90000 - 95 000), very high signal-to-noise {SjN > 500) spectra to determine 
the fraction of odd Ba isotopes (/odd) by measuring subtle asymmetries in the profile of the Ban line at 4554 A. We also use two 
different macroturbulent broadening techniques, Gaussian and radial-tangential, to model the Fe lines of each star, and propagate 
each technique to model macroturbulent broadening in the Ba 4554 A line. We conduct a ID non-LTE (NLTE) treatment of the Fe 
lines in the red giant HD 122563 and the subgiant HD 140283 in an attempt to improve the fitting. We determine [Ba/Eu] ratios for 
the two giants in our study, HD 122563 and HD 88609, which can also be used to determine the relative contribution of the s- and 
r-processes to heavy-element nucleosynthesis, for comparison with /odd. 

Results. We find mathematical solutions of /odd for HD 122563, HD 88609 and HD 84937 of -0. 12 ± 0.07, -0.02 ± 0.09, and -0.05 ± 
0. 1 1 respectively. BD-l-26° 3578 yielded a value for /odd = 0.08 ± 0.08. Only BD-04° 3208 was found to have a physical /odd ratio of 
0. 1 8±0.08. This means that all stars examined here show isotopic fractions more compatible with an s-process dominated composition. 
The [Ba/Eu] ratios in HD 122563 and HD 88609 are found to be -0.20 ±0.15 and -0.47 ±0.15 respectively, which indicate instead 
an r-process signature. We report a better statistical fit to the majority of Fe profiles in each star when employing a radial-tangential 
broadening technique during our ID LTE investigation. 

Conclusions. With the increase of the number of stars for which the Ba isotope fraction /odd has been measured, and the nature of their 
results, there is now a stronger argument to suggest that other synthesis codes that employ alternative approaches to radiative transfer 
(e.g. 3D hydrodynamics) have to be considered to tackle the high level of precision required for the determination of isotopic ratios. 
We have shown that, from a statistical point of view, one must consider using a radial-tangential broadening technique rather than a 
Gaussian one to model Fe line macroturbulences when working in ID. No improvement to Fe line fitting is seen when employing a 
NLTE treatment of the Fe lines. 

Key words. Stars: Population II - Stars: abundances - Galaxy: evolution - Nuclear reactions, nucleosynthesis, abundances 
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Nuclei heavier than the Fe-peak are mainly synthesised via 
two neutron-capture processes, the slow (s-) and rapid (r-) pro- 
cess. For the s-process, neutron-capture rates are much lower 
than y6-decay rates in unstable isotopes, whereas for the r- 
process, neutron-capture rates are higher than y6-decay rates. 
Each n-capture process has a different site for nucleosynthe- 
|1957[ ). Low- to intermediate-mass stars 



sis (Burbidge et al. 



(1M0 <: M < 8M0) evolving along the thermal pulsing 
asymptotic giant branch (TP-AGB) provide the necessary con- 
ditions for the "main" s-process, which is responsible for the 
majority of the s-process elements in the solar-system between 
88 < A < 204 (S neden et al._2008^ , releasing free neutrons via 
'^C(p,y)'^N(jg+Ve)'^C(a,n)'°'o reaction occurring i n the He- 
rich zone in radiative conditions ( |Straniero et al.|1997| . Evidence 



of active s-processing in TP-AGB stars can be seen through ab- 
sorption line detections of short lived ''^Tc and '^Tc isotopes 
( [Smith & Lambert! 1988| and references therein) visible at 4238, 
4262 and 4297 A. 

There are many candidates for r-process sites including 
neutron star explosions (|Imshennik||1992|l, neu tron star surface 
explosions ( Bisnovatyi-Kogan & Chechetkin 1979 1 and neu- 



tron star winds (Panov & Janka 2009 Wanajo et al. 2001 



to name a few. However, presently the most favoured sites for 
the r-process are core-collapse supernovae (Wanajo & Ishimaru 



2006j). Extreme temperatures and run away nuclear processes 
produce an extremely high fluence of free neutrons ( [Wanajo | 
et al.|2003 1, which are necessary for r-process nucleosynthesis. 



The presence of seed nuclei with high n-capture cross- 
sections (o"), such as Fe, is critical for n-capture nucleosynthesis. 
As low- to intermediate-mass stars cannot synthesise nuclides 
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up to the Fe peak, high-cr nuclei must be present at the time of 
the star's formation for the s-process to occur. High-mass su- 
pernova progenitors (M > 8M0) reach high enough tempera- 
tures at the very end of their evolution immediately before the 
supernova explosive phase to synthesise nuclides up to the Fe 
peak. Unlike the s-process, the r-process does not need the pres- 
ence of high-cr nuclei at the time of formation as they are pro- 
duced in situ shortly before the end of the life of the star. Also, 
low- to intermediate-mass stars are long lived in comparison to 
high-mass stars. As such, r-process enrichment from supernovae 
explosions should dominate in the early universe, i.e. in metal- 
poor regimes, with the s-process signatures becoming increas- 
ingly dominant in more metal-rich regimes. This theory was set 
out by Truran| ( 1981[ l and the chemical evolution of n-capture el- 
ements from Ba to Eu was quantified by TravagUo et al. ( 1999 1. 



Indications to support this theory can be seen in evidence pre- 
sented in Frangois et al. (2007, their Fig. 14). 

As the Galaxy becomes more metal-rich over time, s-process 
signatures in stars begin to increase relative to the r-process for 
Ba ( Fran9ois et al. 2007\. W hen one compares the [Ba/Eu] ratios 
from Fran9ois et al. (2007)) with those from Mashonkina e t al.| 
( |2003| l, who study the metallicity at which the s-process begins 
to increase relative to the r-process for the halo and thick disk re- 
spectively, there are variations in star-to-star compositions which 
lead to diff'erent r- and s-process regimes for a given metallicity. 
In particular jFrangois et al. ( |2007 ) find that the s-process begins 
to increase relative to the r-process at [Fe/H] > -2.6, whereas 
[Mashonkina et al.| ( |2003] l find this to occur at [Fe/H] > -1.5. 

One potential way of detecting r- and s-process signatures 
in a star is to measure the isotopic fractions in heavy elements 
using the profile of their absorption lines. There are differences 
between pure s- and r-process isotope ratios in most heavy el- 
ements, which can be in principle detectable through minute 
changes in line asymmetry. Ba is an attractive heavy element for 
which to use of this method, as the hyperfine splitting (hf s) of its 
4554 A line from the singly-ionised stage is quite large ( Rutten 
^ and it offers the possibility of measuring the odd fraction 
via resolved asymmetric lines. 
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The r- and s-process are responsible for five of the seven 
stable isotopes of Ba (the Ughtest two, '^" '^^Ba, arise in the 
so-called p-process (Burbidge et al. 1957\ ). Whereas the s- 
process can synthesise all five Ba n-capture isotopes, shielding 
by '^"^ '^^Xe prevents the r-process from synthesising two of the 
even isotopes, '^'*''^*Ba. 

Using nucleosynthesis calculations ( [Arlandini et al. 1999[ l, 
the values of /odd can be determined for the r- and s-process. For 
a fully s-process regime, /odd.j - 0.11 ± 0.01, and in a fully r- 
process regime, /odd.r = 0.46 ±0.06. Gallagher et al.| ( |20l0 ) show 
a linear relationship between /odd and r- and s-process contribu- 
tions for Ba determined from Arlandini et al. ( 1999| l. Values of 



0.0 < /odd < 0.11 or 0.46 < /odd ^ 1.0 are not physical in the 
context of the nucleosynthesis model but are achievable from 
more ad hoc isotopic mixes. However, we have assumed that the 



theory in |Arlandini et al. 



1999 1 is accurate and we state through- 



out our paper that any value of /odd which Ues outside the limits 
0. 1 1 < /odd ^ 0.46 is non-physical. 

From the point of view of spectroscopy, the even Ba iso- 
topes contribute principally to the formation of the line centre 
in the Ban 4554 A line. The odd isotopes, which are hyperfine 
split, contribute to the spectral region closer to the wings of the 
line (see Fig.[T]i. The relative strength of the odd isotopes located 
toward the blue wing are smaller than those located toward the 



red wing and are further from the line core. As such, when the 
odd isotope contribution to the total line strength is increased, the 
asymmetry in the absorption line's profile is increased. When Ba 
is dominated by the s-process, more of the abundance is associ- 
ated with the even isotopes, which contribute to the line near its 
centre; the line profile has a deeper core with shallower wings 
and the line's asymmetry is reduced. 

There are, however, difficulties in determining isotopic frac- 
tions in Ba, chief among which is the high precision required 
in the analysis of /odd, which demands the highest quality ob- 
servations and radiative transfer models. Issues then arise when 
complex astrophysical behaviours, such as convection, start to 
become visible in the high quality stellar data. Basic assumptions 
used in conventional spectrum synthesis codes that assume a 
plane-parallel geometry (ID) and local thermodynamic equilib- 
rium (LTE) cannot actually replicate the observed Ban 4554 A 
line profile (Galla gher et al.|2010| l. 

The well studied metal-poor subgiant, HD 140283 illustrates 
the difficulties that are encountered when determining /odd, 
which has been attempted several times for this star. Gallagher] 
[eTaL] ( [20T0l l find [Fe/H] = -2.59 ± 0.09 and a low [Ba/Fe] 
ratio - -0.87 + 0.14, which seemingly point to an r-process 
origin as Ba is mainly an s-process element in the solar sys- 
tem. Yet according to their isotopic analysis they find /odd = 
0.02 ± 0.06, which agrees with the result published by Magain] 
(1995) , who finds /odd = 0.08 ± 0.06. Both isotope fractions in- 
dicate a fully s-process regime, contradicting the Truran ( 1981[l 
model. However, measurements of /odd in ID LTE by 'Lambert 
& AUende Prieto ( |2002| l and (for the same spectrum) Collet 



et al. (2009 1 show HD 140283 to be slightly r-process dominated 
with /odd = 0.30 ± 0.21 a r id /odd = ( ).33 + 0.13 respectively , 
which contradicts Magain (19951 and [Gallagher et aL (2010 1. 
When Collet et al. ( 2009!) reanalysed the same spectrum again 
using 3D hydrodynamic model stellar atmospheres they found 
/odd - 0.15 ± 0.12, indicating an s-process signature and sup- 
porting results from Magain ( 1995 1 and Gallagher et al. ( 2010[ l. 
If the three studies showing HD 140283 to have an s-process 
signature are correct, the contradiction with Truran's paradigm 
could be explained by how s- and r-processes vary with [Fe/H], 
or by the inhomogeneity of the interstellar medium (ISM) when 
this halo star formed, in which case star-to-star variations of n- 
capture signatures would be common. The cause of the discrep- 
ancy between the results of the various studies of HD 140283 is 
unclear. 

Unlike Ba, Eu in the solar system has a predominantly r- 
process contribution, which Arlandini et al. ( |1999| calculates 
to be 94% of the total Eu. Both stable Eu isotopes, '^'•'''^Eu, 
show significant r-process contributions relative to the s-process 
and these occur in almost equal amounts at a ratio of 0.48:0.52 
for 151:153 ( Arlandini et al.|[T999| l. High abundances of Eu in 
metal-poor stars indicate a strong r-process signature ( |Spite &| 
Spite 1978 , Sneden et al. 2008 1. Therefore it is common practice 
to use the [Ba/Eu] ratio as an indicator of a star's r- and s-process 
ratio ( .Burris et al.|2000t [Honda et al.|2006l [20071 [Sneden et al.| 
|2008l ). 

Limits on [Ba/Eu] abundance ratios can be set for pure s- 
and r-processes and are found to be +1.45 and -0.81 respec- 
tively ( [Burris et al. |2000|l, which were calculated using solar 
abundances in [Anders & Grevesse[([1989| l. Using the theoretical 
abundances in |Arlandini et al.[ ( |1999[ l, we find the [Ba/Eu] limits 
for the s- and r-processes to be -1-1.13 and -0.69 respectively. 



Gallagher et al.,(2010) highlighted the issues that arise when 



/odd = [A'C^^Ba) + iV('"Ba)]/A'(Ba) 



fitting Fe lines assuming ID LTE, particularly in the wings of 
the line, and plotted the average residual for all the Fe lines anal- 
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ysed. They demonstrated the asymmetries that occur in line for- 
mation; ID LTE radiative transfer codes cannot replicate asym- 
metries. Gallagher et al. (2010 1 speculated that a spectroscopic 
analysis based on 3D radiation-hydrodynamic model stellar at- 
mospheres may resolve these problems. In addition they found 
irregularities between the observed and synthetic profiles in the 
Fe line's core. 

In this paper we determine the isotopic fractions of Ba in a 
further five metal-poor stars under the assumption of ID LTE. 
We describe the observations in ^ and the ID LTE analysis of 
the Ba 4554 A line in ^ Due to the lack of metals in metal-poor 
stars, i.e. stars with [Fe/H] < -2, electron number densities are 
low, which drives down opacities in the atmosphere of late-type 
stars. Therefore assuming LTE in line forming regions of the 
stellar atmosphere is no longer valid ( |Mashonkina et al.|2008| l. In 
^of the work presented here we test whether the irregularities 
seen in the Fe line residuals are due to LTE departures, using 
a non local thermodynamic equilibrium (NLTE) treatment for 
HD 140283 and HD 122563. We discuss the results in ^ 

2. Target selection and observations 

The stellar spectra used in this study are almost the highest qual- 
ity spectra of very metal-poor stars obtained using the High 
Dispersion Spectrograph (HDS) ( jNoguchi et al | 2002| l at the 
Subaru 8.2m Telescope. All have high resolution (R = A/AA - 
90 000 - 95 000, calculated from the widths of several hundred 
ThAr lines), and high signal-to-noise {S /N > 500 per pixel, as 
measured around 4500 A). Such spectra are essential for an ac- 
curate measurement of /odd- Specifics of the observations can be 
found in Table [T] 



( |2007l l and |Mashonkina et al.| ( |2008| ). |Mashonkina et al.| ( |2008l l 
complement the [Ba/Eu] determination nicely in HD 84937 cal- 
culating /odd - 0.43 ± 0.14, indicating an almost fully r-process 
regime. However they found /odd in HD 122563 to be 0.22+0.15, 
a mostly s-process regime, which contradicts the [Ba/Eu] abun- 

( |200S| 



see 



dances found in their study and by Honda et al 
Table |2] We also calculate /odd for two more turn-off stars 
BD-04° 3208 and BD+26° 3578, neither of which have previ 
ous Ba and Eu analyses. 



Table 2. Results from previous ID LTE studies to determine 
[Ba/Eu] and/or /odd for stars studied in this work. 



Star 


[Fe/H] 


[Ba/Eu] 


/odd 


Reference 


HD 140283 


-2.59 


> -0.66 


0.01 + 0.04 


(1) 




-2.40 


~ -1.05 


0.30 + 0.21 


(2) 




-2.50 




0.33 + 0.13 


(3) 




-2.70 




0.08 + 0.06 


(4) 


HD 122563 


-2.77 


-0.50 




(5) 




-2.53 


-0.41 


0.22 + 0.15 


(6) 


HD 84937 


-2.15 


-0.69 


0.43 + 0.14 


(6) 


HD 88609 


-3.07 


-0.48 




(7) 



(1) 


Gallagher et al. 


(2010 


1, as measured by the 4554 A line. (2) 


ILambert & AUende Prieto (2002|i. (3) Collet et al. (2009|l. (4) 


Magain 


( 1995J. (5) 


Honda et al. 


P006 1. (6) 


Mashonkina et al.| 



Table 1. Details of the observations of the stellar spectra. 



3. 1 D LTE analysis 



Star 


Date 


Exp. Time (min) 


S/N 


R 


HD 140283 


22/07/01 


82 


1100 


90000 


HD 88609 


20/04/04 
19/10/05 
20/10/05 


210 


750 


90000 


HD 122563 


30/04/04 


90 


850 


90000 


HD 84937 


22/03/03 


180 


630 


95 000 


BD+26° 3578 


17/05/05 


130 


550 


95 000 


BD-04° 3208 


18/05/05 
19/05/05 


180 


580 


95 000 



Table |2] shows previous results on Ba published for four of 
the stars studied in this paper. It further illustrates the difficulties 
in determining /odd, in that for HD 122563, /odd and [Ba/Eu] do 
not support each othe:[^ while for HD 140283 there are large 
discrepancies between /odd determinations. 

The two giants in our study, HD 122563 and HD 88609, and 
one of the turn-off stars, HD 84937, show a strong indication that 
Ba should be r-process dominated based upon [Ba/Eu] abun- 
dance ratios calculated in [Honda et al.| ( |2006| l, [Honda et al.| 

^ One might take the view that these results are consistent with a 
mixed heavy element origin. 



In this section we briefly review the method used to deter- 
mine the isotopic fractions, abundances and broadening values 
of the stars in our sample. All synthetic spectra in this sec- 
tion were created using the ATLAS ( [Cottrell & Norris[|l978[ l 
radiative transfer code with KURUCZ06 model atmospheres 
(http://kurucz.harvard.edu/grids.htinl ). For a more 
extensive description of the processes involved in the following 
procedure, we refer the reader to jGallagher et al. ( 2010) ). Results 
are provided in Table |3] 

3.1. Barium line lists 

Line lists with all hfs components of the Ba line at 4554 A, 
which is used to determine /odd, were constructed using isotopic 
Arlandini et al. ( 1999|l and hfs information 



information from 



from [Wendt et al.[^84| i and [Villemoes et al.[ ( [1993| ) for pure 
s- and pure r-process mixtures (corresponding to /odd - 0.11 
and 0.46 respectively). Hybrid line lists for -0.24 < /odd < 0.46 
were created from these by adjusting the line strengths of the Ba 
isotopes. Further details on this can be found in Gallagher et al. 
( [2010[ l but we remind the reader that /odd - 0.11 is the lowest 
value of /odd achieved in the s-process of [Arlandini et al. ( 1999| l; 
/odd - 0.00 is the lowest value of /odd achievable in an even-only 
isotope mix, and that values of /odd < 0.00 are non-physical, 
mathematical solutions only. 

Isotopic abundances presented in [Arlandini et al. ( [1999| l 
which we use to deter mine /odd.r ™d /odd.i were normalised (by 
Arlandini et al. 1999 1 to the s-process-only isotope '^'^Sm. We 
explored a renormalisation to '^"*Ba and '^^Ba, but found little to 
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no change in fodd,/- 0.49 (renormalised to '''"^Ba) and 0.46 (renor- 
malised to '■'^Ba). The renormalisation does not affect /odd.i, 
which remains /odd,.? - 0.11. This does not significantly alter 
the interpretation of the results presented here. 

No attempt was made to determine /odd for the 4934 A line 
because, as showed by IGallagher et alj (|2010|), analysis of this 
line is extremely difficult and yields large errorbars due to Fe 
blends found in the wings of this line. Nor do we attempt to 
study higher excitation lines of Ba, as their hyperfine splitting is 
much smaller than that at 4554 A. 



prescription we adopt for ^rt assumes equal speeds in the radial 
and tangential flows, and equal temperatures ( |Gray|2008 1. 

The model spectra were synthesised with ATLAS using a grid 
of values ranging from4.00kms"' to 9.00kms"' in steps of 
0.2 km s ' . The instrumental broadening, determined from ThAr 
lines, and represented as a Gaussian, was also included in the 
synthesis. Each line was fit using the same;^'^ code that was used 
in ^3.2. 1| Results for ^rt can be found in Table |3] (rows (17) to 
(24)), where we have also listed the number of lines that are best 
fit by the ^rt approach for each star. 



3.2. Determination of fhe macroturbulence 



To recap on Gallagher et al. ( 2010|l, w e use ax^ code (derived 



from that of Garcia Perez et a. 



( 2009[ l) to compare the observed 



and synthetic spectra of a number of Fe i and Fe ii lines, com- 
puted for different abundances A(Fe2 wavelength shifts AA, 
and line broadening parameters. Two different broadening ap- 
proaches are used, as described below. We derive, for each Fe 
lines analysed, the best fitting abundance, wavelength shift and 
broadening. An ordinary least squares (OLS) fit is calculated 
through the wavelength-dependence of broadening and A(Fe) 
values to determine the best values at 4554 A. 



3.2.1. Gaussian broadening 

One approach to the macroscopic broadening is to adopt a 
Gaussian of FWHM - v^onv, representing the convolution of 
a Gaussian instrumental profile with a Gaussian macroturbu- 
lent profile. As Lambert & Allende Prieto ( 2002| l show, /odd 
is extremely sensitive to Vconv They found ^fndri/^Vr nnv - 
-0.51 (kms-')-i for HD 140283. [Gallagher et aT] {20101 1 found 
an even larger sensitivity, -0.71 (kms~0~^. The effects of this 
large sensitivity can be reduced by increasing the number of Fe 
lines, A^, used to constrain Vconv, as we take the error in Vconv as 
the standard error, cr/ VA^, where cr is the standard deviation of 
I'conv values. In our work, only Fe lines of comparable equiva- 
lent widths (W) to the Ban line are selected, so they would have 
similar formation depths to the Ba line, implying that macrotur- 
bulent effects on Ba would be well described by Fe. We also do 
not analyse strong lines where uncertain pressure effects in the 
line broadening become significant. 

The adopted atmospheric parameters for each star are listed 
in Table |3] (rows (1) to (5)) and the derived broadening is given 
in row (7). 



3.2.2. Radial-tangential macroturbulent broadening 



In Gallagher et al.| ( 20T0| we compared three different broaden- 
ing techniques to find out which of them best fit the Fe Unes. 
We found that using a rotational broadening mechanism, where 
the macroturbulent broadening of the star was represented by 
vsin/ only (with the instrumental broadening still represented 
by a Gaussian), almost always gave a worse fit than when we 
employed a simple Gaussian. However we found that using a 
radial-tangential macroturbulent broadening technique (^rt) al- 
lowed us to fit spectral lines slightly better than the Gaussian 
mechanism. As a result, we have continued to employ this broad- 
ening type in the current investigation for all five stars and again 
for HD 140283, as well as the simple Gaussian approach. The 



3.3. Tile isotopic fraction of barium 

Once Vconv and ^rt were obtained for the two broadening mech- 
anisms, the Ban 4554 A line was analysed using ax^ code very 
similar to the one described in 5 3.2 to find values for the free pa- 
rameters AA, A(Ba) and /odd that minimise for the 4554 A line 
in each star. Results for /odd for Gaussian and radial-tangential 
broadening techniques can be found in Table |3] The best fit Ba 
profiles are shown in Fig. [T] (Gaussian) and in Fig. |2] (radial- 
tangential). From the Gaussian results we find for HD 122563, 
HD 88609, HD 84937, BD-04° 3208 and BD-i-26° 3578 that 
/odd = -0.12 ± 0.07, -0.02 + 0.09, -0.05 + 0.11, 0.18 + 0.08 
and 0.08 + 0.08 respectively. These results would suggest that all 
stars examined here show a high s-process fraction. 



3.4. Tlie [Ba/Eu] ratio 

In our sample only the two giants had Eu line strengths ade- 
quate to conduct a Eu abundance analysis. This was done by 
examining the Euii 4129 A line. The Eun 4205 A line is not 
used in this study as it is blended with a V ii line ( [Honda et al.| 
2006 Gallagher et al.|2()T0 1 that can affect abundance determi- 
nations. Also no attempt is made to analyse the isotopic splitting 
of '^'"'^^^Eu; we assume a fixed 50:50 isotopic split of 151:153 
when constructing the hfs-affected Eu 4129 A line list. The Eu 
line list for the solar r- and s-process ratio was constructed using 
hfs information from |Becker et al.[ ( [1993|l and[Krebs & Winkler[ 
(1960) with \og gf values taken from Biemont et al. ( 1982[ l. 
Abundances for Eu are taken as those that satisfy the x"^ min- 
imum. From these abundances, and those found by the Ba anal- 
ysis, [Ba/Eu] was calculated and the results for both stars can be 
found in Table [3] row (10). 



3.5. Error analysis 



In Gallagher et al.[ (!2010) we demonstrated how altering one stel- 



to 



lar parameter can force other parameters, in particular 
compensate for its effect on /odd. In that paper we called this 
case 1. The compensation by Vconv was well illustrated when one 
looked at the effects of changing the microturbulence (^). Once 
Vconv was recalculated from the Fe lines, the effect ^ had on /odd 
was nullified. As Vconv partially compensates for other changes 
in the stellar parameters as well, i.e. Teff and logg, their effect 
on /odd is also reduced. 

To calculate the error in /odd we look at five possible sources 
of error: Vconv, ^, T^eff, logg and the Unsold approximation en- 
hancement factor, Ey, which enhances the effect of ye.vdw in the 



Van Der Waals calculation. 



76 



T6,vdw£^y. In our analysis we 



have used Ey = 2.2. To test the effect it has on Vconv and /odd we 
have decreased this to 1.5. The effect of uncertainties in [Fe/H] 



A(X): 



log.o(^)+12. 



on /odd is negligible and as such was ignored (Gallagher et al. 
[2010[ ). It is expected that every star belonging to the same stage 
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HD1 22563 



HD88609 




4554.0 4554.1 
\ (A) 



4554.0 4554.- 
\ (A) 



Fig. 1. The best fit Ban 4554 A lines for each star, using a Gaussian broadening technique. Each figure displays the observed Ba 
profile (diamonds) and the best fit synthetic profile (solid line), which includes the error on /odd (dashed line). We have also included 
a schematic of the odd and even isotopes for reference. The lower panel of each figure shows the residuals (obs-syn) of each fit as a 
percentage. For reference we have included HD 140283, which was analysed in [Gallagher et aL] ( |2010| l. 
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HD' 22563 



HD88609 




4554.2 4554.3 



4554.0 4554.1 
X (A) 



4554.0 4554.1 
\ (A) 



4554.2 4554.3 



Fig. 2. The best fit Ba ii 4554 A lines for each star, using a radial-tangential broadening technique. Each figure displays the observed 
Ba profile (diamonds) and the best fit synthetic profile (soUd line), which includes the error on /odd (dashed line). We have also 
included a schematic of the odd and even isotopes for reference. The lower panel of each figure shows the residuals (obs-syn) of 
each fit as a percentage. 
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Table 3. Results of the ID LTE analysis conducted on five metal-poor stars and the various parameters used in their analysis. 





Parameter 


HD 122563 


HD 88609 


HD 84937 


BD-04° 3208 


BD-h26° 3578 


HD 140283 


Atmospheric Parameters 


(1) 


T'eff (K) 


4570 ± 100 


4550 + 100 


6290 + 100 


6340 ± 100 


6240 + 100 


5750 + 100 


(2) 


log g (cms-') 


1.1+0.3 


1.1+0.3 


3.9 + 0.3 


4.0 + 0.2 


3.9 + 0.2 


3.7 + 0.1 


(3) 


[Fe/H] 


-2.77 + 0.19 


-3.07 + 0.20 


-2.15 + 0.30 


-2.28 + 0.20 


-2.33 +0.20 


-2.50 + 0.20 


(4') 


(km s^h 


2.2 + 0.5 


2.4 + 0.5 


1.2 + 0.3 


1.5 + 0.2 


1.5 + 0.2 


1.4 + 0.1 


(5) 


Atmosphere Reference 


(1) 


(2) 


(3) 


(4) 


(4) 


(5) 








Gaussian 


broadening 








(6) 


J oaa 


-0.12 ± 0.07 


-0.02 + 0.09 


-0.05+0.11 


0.18 + 0.08 


0.08 + 0.08 


0.02 + 0.06 


(7) 


Vconv (km s"') 


6.99 + 0.07 


7.03 + 0.08 


6.98 + 0.06 


6.62 + 0.05 


6.41 + 0.06 


5.75 + 0.02 


(8) 


[Fe/H] 


-2.90 ±0.17 


-3.17 + 0.17 


-2.24 + 0.11 


-2.42 + 0.11 


-2.45 + 0.12 


-2.59 + 0.09 


(9) 


[Ba/Fe] 


-0.85 + 0.21 


-0.91+0.21 


0.04 + 0.15 


-0.20 + 0.15 


-0.06 + 0.16 


-0.87 + 0.14 


(10) 


[Ba/Eu] 


-0.20 + 0.15 


-0.47 + 0.15 








> -0.66 


(11) 


<A^) (mA) 


-15.86 + 0.37 


-14.70 + 0.90 


-14.69 + 0.53 


-11.47 + 1.093 


-13.29 + 0.89 


-12.67 + 0.48 


(12) 


5/odd/^Vconv (kms"')"' 


-0.57 


-0.74 


-0.89 


-0.78 


-0.77 


-0.71 


(13) 


5/oddM^(kms-i)-i 


-0.08 


-0.08 


0.05 


0.05 


0.05 


0.00 


(14) 


5/oddMreff (100 K)-i 


-0.04 


-0.05 


-0.03 


-0.02 


-0.02 


-0.02 


(15) 


<5/odd/51ogg (dex)"' 


0.02 


0.02 


-0.27 


-0.23 


-0.23 


-0.28 


(16) 


Sfodd/SEy (0.7)-i 


-0.02 


-0.02 


-0.05 


-0.04 


-0.04 




Radial Tangential broadening 


(17) 


/odd 


-0.16 + 0.05 


-0.11 +0.07 


-0.02 + 0.10 


0.20 + 0.10 


0.09 + 0.10 


0.02 + 0.04 


(18) 


^RT (kms"') 


5.69 + 0.09 


5.74 + 0.09 


5.77 + 0.08 


5.36 + 0.06 


5.14 + 0.07 


4.30 + 0.02 


(19) 


Vinst (kms"') 


3.32 


3.32 


3.20 


3.20 


3.20 


3.32 


(20) 


number of best fit Fe lines 


53 


34 


30 


28 


32 


38 


(21) 


^/oddM^RT (kms"')"' 


-0.30 


-0.64 


-0.60 


-0.73 


-0.70 


-0.63 


(22) 


5/oddM^(kms-i)-i 


-0.08 


-0.08 


0.05 


0.05 


0.05 


0.00 


(23) 


^/oddM^eff (100 K)-i 


-0.01 


-0.001 


-0.05 


-0.07 


-0.06 


0.01 


(24) 


5/odd/51ogg (dex)"' 


-0.01 


-0.02 


-0.18 


-0.22 


-0.21 


-0.35 


Other spectral information 


(25) 


Spectral range (A) 


3080 - 4780 


3070 - 4780 


4130 - 6860 


4130-5340 


4050 - 5250 


4100 - 6900 


(26) 


Fe line sample size 


54 


35 


44 


36 


37 


93 


(27) 


WBa (mA) 


99.1 


93.2 


49.7 


34.7 


42.7 


20.1 



(l) |Honda et al] ( [20061 ) . (2) [Honda et aL] ( [2007) . (3) |Aokiet al.] ( [2009[ ). (4) [Garcfa Perez et al.| ( [2009[ ). (5) [Gallagher etaL[ ( [20T0j i 



of evolution, i.e. giant, sub-giant and turn-off, would reproduce 
comparable sensitivities to each stellar parameter. As such we 
have run each test for two of the five stars: to test the sensitivity 
of /odd in the giants we have used HD 88609, and for the turn- 
off stars we have chosen BD-04° 3208. However, we have run 
sensitivity tests of Vconv and ^rt for all stars as these parame- 
ters have the largest effect on /odd- The results of the sensitivity 
of /odd can be seen in Table [3] rows (12) to (16) and rows (21) 
to (24). The tests confirm that the /odd determinations are es- 
sentially unaffected by the choice of Ey, Teff or ^, but log g and 
the macroturbulent broadening effects are more important. It has 



been reported by |Tajitsu et al.[p010j ) that the EEV 42-80 CCDs 
used in the HDS suffer from nonlinearity. By investigating the 
differences between the corrected and non-corrected spectra, we 
found that its effect on /odd was negligible. 



4. 1D NLTE Fe line analysis 



In Gallagher et al. ( 2010[ ) we aimed to constrain the macroturbu- 
lence of HD 140283 through the use of the Fe lines. Using the 
same method, we have also done this for five further stars in ^ 



This was done under the assumptions of LTE. However, it is well 
known that Fe i suffers from the effects of NLTE (Thevenin &| 
[Idiart 1999 Shchukina et al. 2005 i in metal-poor stars. We there- 
fore sought to quantify NLTE effects for Fe on the preceding 
analysis, in particular on the determination of macroturbulence, 
and therefore the value of /odd- No attempt is made to determine 
NLTE corrections for Ba itself, either in /odd or the [Ba/Eu] ra- 
tio. [Mashonkina et al.] ( [2008| l demonstrated that corrections to 
the abundance ratio are not significant enough to change the in- 
ferred r- and s-process regime in a star. The determination of /odd 
assuming NLTE goes beyond the scope of the work presented in 
this section. 

To compute the Fei profiles, a version of the NLTE code 
MULTI ( |Carlsson]|198 6]l was used with modifications to include 
the effects of line-blanketing, described in [Collet et al. (2005Jl. 
The code employs MARCS model atmospheres. The model atom 
used was that adopted by Hosford et al. ( 2010| l. For a longer dis- 
cussion on the model atom, processes of NLTE radiative transfer 
and its impact on stars of this type, see Hosford et al.[ ( 20T0| . It 
is important, however, to mention the parameter Sh, the scahng 
factor for the collisions due to H as described by the approxi- 
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mate Drawin formula ( |Drawin|[l968| |1969| l. Due to the uncer- 
tainties in the magnitude of the H collisions, the value of 5h 
is still uncertain and is treated differently by different works. 
[Collet et al.| ( |2005 i treats it as a free parameter and tests val 



Table 5. Results from the NLTE analysis and the ATLAS LTE 

results for comparison. 



ues 5h = 1 and 0.001. Korn et al. ( 2003[) found a higher value, 
5h = 3, however in more recent work ( |Mashonkina et al.|20T0| l 
the same group has constrained it to 0.1 based on an improved 
Fe atom. Here we adopt the values 1 and 0.001; this gives us 
two sets of synthetic spectra, one with the Drawin (Sh - 1) de- 
scription of H collisions, and a second close to maximal NLTE 
effects (5h = 0.001). We also compute line profiles in LTE using 

MULTI. 

NLTE calculations of Fe i line profiles were performed for 
two stars, the subgiant HD 140283 and the giant HD 122563. 
The first was analysed so a comparison could be made between 
this work and the LTE analysis of Gallag her et al.| ( |2010[ ). The 
second was analysed to understand the NLTE effects on Fe i in 
giants, where the atmospheres are more tenuous and also cooler, 
to act as a comparison to the LTE work in this paper. 

Three sets of MULTI runs where completed for each star, one 
in LTE, and two in NLTE with the S h values mentioned above. 
A analysis procedure similar to the one described in ^been 
used. Synthetic profiles were created over the parameter space 
shown in Table |4] The extrinsic broadening was represented by 
a Gaussian consolidating both macroturbulence and instrumen- 
tal broadening, the values FWHMg in Table |4] representing the 
FWHM of the Gaussian. 



Table 4. Ranges of parameters used to create the synthetic pro- 
files in MULTI. 







grid parameter ranj 


;es 






A(Fe) 


FWHMg 




Star 


Run 




(mA) 


(mA) 


HD 140283 


LTEmulti 


4.28-5.28 


73-119 


+25 




Sh^I 


4.66 - 5.28 


73-119 


+25 




Sh = 0.001 


5.00-5.70 


73-119 


+25 


HD 122563 


LTEmulti 


4.00 - 5.58 


74-114 


±25 




Sh^I 


4.00 - 5.58 


74-114 


+25 




Sh = 0.001 


4.68 - 5.50 


74-114 


+25 



The intention was to see if the observed Fe i profiles could be 
better fit by spectra computed in NLTE, in particular for the gi- 
ants, where the cooler atmospheres result in larger line strengths 
(at fixed [Fe/H]). Here, the synthetic LTE core is too shallow and 
the wings too broad even at the minimum value of the line. In 
switching to NLTE, we hoped we might model the Fe lines bet- 
ter, and hence obtain better values of macroturbulent broadening 
and, in later works, a more reliable value of /odd for Ba. 

Table [s] gives the results from the x^ analysis using MULTI, 
and the atlas LTE results for comparison. The MULTI values are 
the mean of 51 lines for HD 140283 and 31 lines for HD 122563; 
errors are the standard error. Fewer lines are used in the MULTI 
analysis than in the ATLAS LTE analysis due to the incomplete- 
ness of the adopted model atom which leads to several of the 
lines not being computable in MULTI. Also, after further scruti- 
nizing the results, three lines from HD 140283 and three lines 
from HD 122563 were removed from the analysis. This was due 
to either a poorly defined continuum, or blending lines very close 
to the line of study. 



Star 


Run 


<A(Fe)) 


(VconvXkmS ') 


HD 140283 


LTEmulti 


4.96 + 0.01 


6.03 ± 0.04 




Sh^I 


5.28+0.01 


5.98 ± 0.05 




Sh = 0.001 


5.53 + 0.01 


5.92 ± 0.04 




LTE ATLAS 


4.92 + 0.04 


5.76 ± 0.09 




(51 line subset) 








LTEatlas 


4.91+0.01 


5.75 ± 0.02 




(93 lines) 






HD 122563 


LTEmulti 


4.58 + 0.03 


7.02 + 0.10 




5h = 1 


4.98 + 0.03 


6.95 ± 0.09 




Sh = 0.001 


5.12 + 0.03 


6.91+0.09 




LTEatlas 


4.54 + 0.11 


7.06 + 0.18 




(31 line subset) 








LTEatlas 


4.51+0.04 


6.99 ± 0.07 




(54 lines) 







A comparison of MULTI LTE and ATLAS LTE results shows 
that the abundances and extrinsic broadening value are compa- 
rable for both cases, the only exception being the broadening 
values of HD 140283, which is the most important parameter as 
it directly affects /odd. For HD 140283 we have larger values in 
NLTE than that of [Gallagher et al.| ( |20T0| . Table |3] shows that 
for HD 140283, choosing the same 51 line subset doesn't "fix" 
the discrepancy. For HD 122563 the difference is much smaller 
and is within the errors. One could naively ascribe this differ- 
ence to be discrepancies in the MARCS and KURUCZ model atmo- 
spheres, e.g. different temperature gradients and density gradi- 
ents. However, a comparison between the MARCS and KURUCZ 
model atmospheres can be found in| Hosford (2010 ), who found 
little to no effect by these parameters. It is most likely that the 
differences arise from the different treatments of the broadening 
itself in MULTI and ATLAS. We see comparable wavelength shifts 
between the two analyses, and comparable profiles for each line. 
In NLTE we see an increase of abundance with decreasing Sh 
due to the overionisation effects becoming more pertinent; thus 
a larger positive abundance is needed for neutral lines to com- 
pensate. As a consistency check, NLTE abundances from the;if^ 
analysis were compared to those from an equivalent width anal- 
ysis and were found to be comparable at each Sh value. 

An interesting thing to note here is the sensitivity of Vconv to 
the value and the difference between LTE and NLTE values. 
For both stars, Vconv decreases by 0. 1 km s" ' in going from MULTI 
LTE to MULTI^H - 0.001. Naively, one might infer from Table|3] 
that this would increase /odd from an LTE analysis of Ba 4554 A 
between 0.06 and 0.09. However, some of this difference may 
not fully apply to Ba calculated in LTE, due to different broad- 
ening under NLTE. Nonetheless, it points to the importance of 
an accurate description of the radiative transfer if one is to be 
able to determine /odd. 

One hope of the NLTE analysis was that it may aid in fitting 
the stronger Fe lines of the giant by producing more realistic line 
profiles, and hence reduce the residuals in the Fe analysis (see 
Gall agher et al.|20"T0 their Fig. 10). In Fig.[3]we regenerate the 
plot of average residuals for HD 140283 using a MULTI analy- 
sis of 51 Fe lines. We have co-added the residuals from each Fe 
line to smooth out any unique line defects. The residuals are of 
similar magnitude to those from |Gallagher et aL] ( [2010| l and have 
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HD140283 HD122563 




a A (A 



Fig. 3. 

MULTI 



Plots of average residuals from;\f^ analysis of Fe lines in HD 140283 and HD 122563, for Sh - 1, 0.001 and LTE using 
For comparison we have generated the same plot using atlas output data for the same subsets of hnes. 



similar features. Most importantly, the MULTI residuals for NLTE 
with different Sh values and for LTE are all very similar. There 
is in fact a very slight worsening of the core fit with increased 
NLTE effects. It is also seen that the residuals for HD 122563 are 
of a similar form, albeit with a greater magnitude. The lines anal- 
ysed in HD 122563 are far stronger than those in HD 140283. As 
the lines become stronger, they become harder to fit precisely, 
but using NLTE profiles has no effect on improving this. 

The fact that the MULTI NLTE analysis has led to no improve- 
ment over the MULTI LTE analysis in the fit of the Fe lines is not 
entirely surprising. In fact when compared to the ATLAS analysis 
for the same subset of Fe lines (seen in red) there is a consid- 
erable worsening of the residuals, particularly for HD 140283, 
although NLTE effects were still important to investigate. The 
effects of NLTE on Fe i in metal-poor stars are dominated by 
overionisation. To a first approximation, the populations of neu- 
tral Fe energy levels can be seen as following a Boltzmann dis- 
tribution relative to one another, but with a shift in the ionization 
equilibrium. The results show that this does not affect the shape 
of the profiles, though for a given equivalent width the abun- 
dance needed to reproduce the line increases. The other main 
NLTE effect, a change in the line and continuum opacities with 
height in the atmospheres, apparently has little overall effect on 
the Fe i line profiles. 

In summary, NLTE effects in Fe affect Vconv by up to ~ 
0.1 kms"', and may affect /odd for Ba up to 0.06 or 0.09, but 
NLTE does not result in a better fit to the Fe lines (Fig. [3]l. 
Other mechanisms may come into play, e.g. as represented by 
3D model atmospheres. It would still be interesting to see, how- 
ever, how a NLTE treatment of Ba affects the inferred /odd- 



5. Discussion 

With the exception of one star, BD-04° 3208, all of the stars 
analysed in this paper and HD 140283, which was studied in de- 
tail in |Gallagher et al.| ( |2010 ), show a non-physical isotope ratio 
(/odd < 0.11) close to the s-process-only composition. The non- 
physical results for /odd suggest that applying a ID LTE treat- 
ment to analyse the isotopic fraction of the Ba 4554 A line does 
not appear to be very robust. There are several possibilities why 
this might be the case, and a few possible solutions to the prob- 
lem, which we now discuss. 



5.1. Alternative broadening tecliniques 

We have shown that /odd has a high sensitivity to Vconv, which 
we have determined by fitting synthetic ID LTE (and in ^ 
NLTE) profiles to Fe lines. There are several reasons why we 
chose to use Fe lines. Firstly they are the most abundant species 
in a metal-poor spectrum and cover a considerable range in line 
strength. As mentioned above, these two attributes are useful be- 
cause they allow us to use lines that form at similar depths to the 
Ba line. However, we revisit this argument and consider whether 
Fe is a sensible atomic species to use to determine the macrotur- 
bulence of another species. 

To test this hypothesis, we analysed several Ca i lines in the 
star with the best quality stellar spectrum, HD 140283 {S /N = 
1100), using the same techniques described in ^ Only seven 
lines were found to be unblended and to have similar strengths 
to the Ba Une (10 < W (mA) < 50). This meant that the stan- 
dard error in v^onv was much larger for this set of lines than was 
found for the 93 Fe lines we used in Gallagher et al. ( 2010| l. 
We found that (Vconv)ca = 5.63 + 0.13 kms"'. This is a smaller 
value than was determined using Fe lines (5.75 ± 0.02 kms"') 
but still within the 1 cr eiTor, and the statistics of using only seven 
lines meant that by using Ca we would greatly increase the un- 
certainty in /odd- Nevertheless it was found that /odd, which from 
the Fe analysis was found to be 0.01 + 0.06 (in the 4554 A line|^, 
would move higher to 0.10 + 0.1 1. We can not say whether Ca 
is a better atomic species to use than Fe as the spread in Vconv is 
too high. 

To avoid using other elements to constrain the macroturbu- 
lent broadening of Ba, we experimented by treating Vconv as a 
free parameter whilst determining /odd, thus deriving it from the 
Ba line. This was done for all six stars in this study. Results are 
shown in Table|6] In Table|6]we also give the standard deviation 
(s.d.) of the Fe line measurements, as an indication of the un- 
certainty associated with the measurement of just one line, e.g. 
Ba 4554 A. As shown, Vconv is little changed (within 1 s.d.) in 
four stars when using just the Ba line, and within 2 s.d. in a fifth 
star, and therefore there is little change to /odd- It is interesting to 
note that two of the turn-off stars and the subgiant, HD 140283, 
are found to have higher Vconv values, driving smaller /odd (more 
non-physical) ratios though in all cases the differences in Vconv 
are comparable to the uncertainty expected for one measure- 
ment. It is therefore still unclear whether the Fe analysis does 

Table [3] lists f^a as calculated by both the 4554 and 4934 A lines 
for HD 140283. 
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Table 6. Values of Vconv (measured in km s"') as determined from 
Fe lines and directly by the Ban 4554 A line. Note that we can- 
not calculate /odd in BD-i-26° 3578 when we employ Vconv (Ba). 



Star 


^conv 


s.d. 


/odd 


^conv 


/odd 




(Fe) 


(Fe) 


(Fe) 


(Ba) 


(Ba) 


HD 122563 


6.99 


0.54 


-0.12 


6.82 


-0.04 


HD 88609 


7.03 


0.45 


-0.02 


6.98 


0.01 


HD 84937 


6.98 


0.41 


-0.05 


6.92 


-0.03 


BD-04° 3208 


6.62 


0.32 


0.18 


6.88 


-0.01 


BD+26° 3578 


6.41 


0.34 


0.02 


8.85 




HD 140283 


5.75 


0.19 


0.01 


6.06 


-0.23 



describe the Doppler broadening for Ba well or not, but there is 
not strong evidence against it. 

As well as testing whether Fe was an adequate species to use, 
we also considered whether a simple Gaussian adequately de- 
scribes macroturbulent broadening. In Gallagher et al. (2010) we 
showed that a ^rt macroturbulent broadening mechanism better 
fit several Fe lines than a Gaussian. Therefore we employed a 



^RT treatment for each star analysed in this investigation (§3.2.2 
Fig.|2] and Table[3]l. It was found that for the giant stars, a radial- 
tangential fit the Fe lines better, with only one line in each spec- 
trum better fit by a Gaussian. We found that for the three turn- 
off stars a ^rt technique fit 68% - 86% of Fe lines better. This 
demonstrates a clear reason to consider using ^rx when work- 
ing in ID LTE over a Gaussian broadening mechanism. This is 
particularly clear with the giant stars, whose Fe line cores are 
closer to saturation causing the wings to become more signifi- 
cant. Examination of the residual plots in Fig. |4] shows further 
the improvement to the fits when using ^rj. However we stress 
to the reader that whilst using such a technique under the as- 
sumptions of ID LTE seems to improve upon fitting errors asso- 
ciated with using a traditional Gaussian, both are still symmetric 
profiles and are unable to remedy the issue of asymmetries asso- 
ciated with absorption lines in real stellar spectra. 



5.2. Fe line residuals 

Fig. |4] shows the average residuals of the Fe lines used in deter- 
mining Vconv and ^RT. It illustrates the difficulties in fitting ab- 
sorption lines with ID LTE synthetic profiles. 

It is seen that the turn-off and subgiant stars have quite 
an asymmetric residual profile, with particular problems occur- 
ring in the red wings, 100 to 130 mA from line centre for all 
four stars. This seems to be caused by underlying assumptions 
adopted in ID LTE radiative transfer codes; real absorption lines 
are not perfectly symmetric. 

This investigation is an extension of the asymmetry analy- 

'(20TO'); 



sis conducted in Gallagher et al. (2010); the asymmetry seems 
to occur in all four stars. Severe fitting issues occur in the giant 
stars, HD 122563 and HD 88609, where the Ba (and hence Fe) 
line equivalent widths are ~ 90 mA (Tabled row (27)) and lines 
cores begin to saturate so pressure broadening becomes more 
significant in the wings, and this may explain the symmetric 
residuals seen at +0.14 A. 



Gaus&ion 
Radial Tongentio 





Rodial Tongentio 




Radial Tongentio 





Fig. 4. The average residuals of the Fe lines used to constrain 
Vconv and ^RT. 
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5.3. New approaches to determine isotope mixtures 

It is clear that there is still uncertainty over which atomic species 
to use to determine Vconv However it is more likely that the as- 
sumptions that a ID LTE code adopts, in particular the sym- 
metric broadening mechanisms used in replicating macroscopic 
turbulence in a star's atmosphere, is a more likely reason for the 
large residuals, illustrated in Fig. |4] In the previous paper we 
asked, but did not answer, whether 3D hydrodynamical codes 
could solve these issues. One major drawback of using 3D tech- 
niques is the time-scales involved with the both the model atmo- 
sphere production, which cuiTently can take several months of 
computation, and the radiative transfer calculations, which cur- 
rently take several hours per spectral line synthesis. In addition, 
3D hydrodynamics is a rather new field of analysis, and as such 
there has been very little published to test the reliability of the 
codes currently available. However it appears that 3D hydrody- 
namics is the next logical step to take in this field ( [Collet et al.| 
2009 1 Ludwig & Kucinskas 2005 ) in response to the asymmetric 
profiles. However, this may not solve the problem of /odd since 
Collet et al. ( 2009 1 found that while the 3D model atmosphere 
did seem to improve upon asymmetries associated with the Ba 
4554 A absorption line, the result for /odd increased the s-process 
contribution relative to the ID LTE result. Therefore, even if 
adopting a spectral line synthesis code that employs 3D hydro- 
dynamic model stellar atmospheres is the next step to better un- 
derstanding the problems with fitting absorption lines, it does 
not yet appear to be able to reduce the high s-process fractions 
we have found, in fact it would appear using such a technique 
would enhance the s-process fraction further. 

In addition to new approaches in the analysis of isotopic frac- 
tions, perhaps a new, more accurate observing technique is re- 
quired to get the best possible data for this type of complex anal- 
ysis. Astronomers using precise radial velocities to detect extra- 
solar planets via the parent star Doppler "wobble" are employing 
observing procedures to attain Doppler precisions < 3 ms"' by 
passing the starlight through an iodine (I2) absorption cell. This 
allows them to measure very small wavelength shifts, which af- 
fect the instrumental profile ( Butler et al. 1996 1. This can be cal- 
culated over all wavelength scales and in each short exposure, to 
determine the instrumental broadening profile at any time. 

(^04) have shown that there is a very slight 



to-noise (S /N > 550) spectra. We found that all stars show a 
high s-process signature, and only BD-04° 3208 had an /odd 
value (0.18 ± 0.08) that was physical, 0.11 > /odd > 0.46 ac- 
cording to the isotopic abundance determinations by Arlandim] 



et al. 



1999[ l. According to these limits, all other stars yield a 
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asymmetry in the ThAr lines recorded with HDS, at least in the 
HD 140283 observing run, which were used to calculate Vinst. 
Using the I2 method described above, it is possible to decon- 
volve the instrumental profile from the stellar spectrum, leaving 
only real stellar broadening in the lines. This would remove most 
of the uncertainty surrounding the instrumental profile from any 
subsequent analysis. Issues with this method are the requirement 
for a stable I2 cell and the large amount of time that would be 
required to complete such a complex reduction. In conjunction 
with a 3D hydrodynamical treatment this might require a lot of 
time and still may not explain the r- and s-process contributions 
any better than the method described in this paper. However, it 
is a method still favoured by astronomers looking for radial ve- 
locity fluctuations in stars, as it yields the smallest error in the 
spectral data. The difficulty in modelling the Fe lines, let alone 
Ba, may justify the additional effort required. 



6. Conclusions 

We have carried out a careful examination of the Ba isotopic 
ratios for five metal-poor stars using a ID LTE treatment in con- 
junction with high resolution (R = A/AA > 90000), high signal- 



non-physical isotopic fraction, /odd < 0.11. 

Using the radiative transfer code MULTI with MARCS ID 
model atmospheres, we have also applied a ID NLTE treatment 
to Fe fines in HD 122563 and HD 140283. We found tfiat us- 
ing NLTE did not fit the Fe line cores better, rather it increased 
the core residual in both stars, see Fig. [3] A large NLTE effect 
{Sh - 0.001) did, however, reduce Vconv by 0.1 kms"', and may 
imply an increase of /odd by ~ 0.06 - 0.09, which would reduce 
the non-physicality of the results. 

During the investigation we have shown how ID LTE radia- 
tive transfer codes and model atmospheres seem to be inadequate 
to determine /odd- In particular we have shown that asymmetries 
in the red wings are observed in the Fe line residuals for all 
four turn-off' and subgiant stars, while the two giants have larger, 
albeit symmetric, residuals. We speculate that spectral line syn- 
thesis with 3D hydrodynamic model stellar atmospheres would 
improve residuals in both the Fe and Ba lines, but we make 
no predictions as to whether it would improve upon the obvi- 
ous problems with non-physical isotopic fractions we have en- 
countered here. In addition we have asked, but not answered, 
whether taking the observations using an I2 cell, designed to 
improve accuracy and quality of the data, would help to better 
constrain /odd- Certainly if one compares the error estimation of 
/odd in HD 140283 (+0.06, S /N = 1 100) with HD 84937 (+0.1 1, 
S /N - 630), there is almost double the uncertainty in HD 84937 
than in HD 140283. We speculate that this could be due to the 
quality of the spectral data used, however, we acknowledge that 
this could also be a result of the difference in the number of Fe 
lines used to determine Vconv^ as we measure error by the stan- 
dard eiTor. However, crv^„, ,hd 84937 = 0.06 km s"' is fairly typical 
as an error estimate for these stars, yet o-f^^^ is larger in HD 84937 
than in all other stars. 

It was found that using a radial-tangential broadening tech- 
nique rather than a traditional Gaussian improved upon fitting 
errors for both the Fe and Ba lines. Whilst using such a tech- 
nique does not remove asymmetries seen in absorption lines for 
real stellar data, it appears that if one were to adopt a ID LTE 
approach to resolve an absorption line, a radial-tangential profile 
should be used to model the macroturbulence of the star rather 
than a Gaussian, although a Gaussian could still be used to model 
the instramental profile. 

We have also conducted a study of Eu abundances in the 
two giants, HD 122563 and HD 88609. This was used to deter- 
mine [Ba/Eu] ratios, which are used to assess the star's s- and 
r-process content. Both stars have [Ba/Eu] ratios that indicate a 
large r-process contribution, which agrees well with Honda et al. 
(2006) & Hon da et al.| ([20071), but contradicts the isotope anal- 
ysis conducted in this paper. This further justifies scepticism of 
the ID LTE techniques employed in our Ba isotope analysis. 

We have increased the number of stars that have undergone a 
Ba isotope analysis. However it is difficult to believe that all stars 
analysed have a significant s-process contribution, or s-process 
enhancement in contradiction of the (limited) [Ba/Eu] data and 
the strengths of the |Truran[ ( [1981| theory and jTravaglio et al.[ 
(1999) calculations. It is much more likely that the symmetric 
ID LTE techniques used in this investigation are inadequate and 
improvements to isotopic ratio analysis need to be made. 
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